Supershear Rayleigh waves at a soft interface 
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We report on the experimental observation of waves at a liquid foam surface propagating faster than the bulk 
shear waves. The existence of such waves has long been debated, but the recent observation of supershear 
events in a geophysical context has inspired us to search for their existence in a model viscoelastic system. An 
optimized fast profilometry technique allowed us to observe on a liquid foam surface the waves triggered by 
the impact of a projectile. At high impact velocity, we show that the expected subshear Rayleigh waves are 
accompanied by faster surface waves that can be identified as supershear Rayleigh waves. 



Supershear rupture, i.e., rupture propagating at unusual 
speed, faster than that of shear waves, has recently been ob- 
served in large seismic events LIJ and experiments |2|. These 
events break the theoretical limit stating cracks should not 
propagate faster than the velocity of Rayleigh waves (trans- 
verse surface waves) |3 |. Their increased hazard raised 'the 
need to study speed' B | and they were successfully explained 
using both internal friction and soils heterogeneities. Inspired 
by these remarkable observations, we aim to address a related 
question: could Rayleigh waves 1 5 1 themselves propagate at a 
velocity faster than the shear waves? The existence of super- 
shear Rayleigh waves has long been debated, either in pure 
elastic I5H3 or dissipative media |[Tq1IT21| . but lacks experi- 
mental observations. 

The model viscoelastic material that we use to address ex- 
perimentally the aforementioned question is a liquid foam 
1131 . It consists of a concentrated dispersion of gas bubbles in 
a continuous liquid phase. With its low shear modulus, caus- 
ing waves to propagate slowly, foam is ideal for wave velocity 
measurements. Furthermore, because of its low density, iner- 
tia is negligible. However, strong dissipation makes acoustics 
methods challenging in viscoelastic media such as foams lfT4l . 
and, although the rheology of these materials is widely stud- 
ied fTJ', 15], little is known about the propagation of surface 
waves in foams. Using the impact of a solid bead, we trigger 
deformations of a foam surface, that we measure by a fringe 
projection profilometry technique L16, 17 1 already tested for 
fast impacts on liquids ifTSl as shown in Fig. [TJa,b). We 
demonstrate that this non-invasive technique can accurately 
track impact-generated surface waves. For gentle impacts, our 
measurements are in agreement with the speed predicted from 
published rheology experiments 1 19 |. But the foam response 
changes as impact velocity increases, and some surface waves 
are found to propagate faster than shear waves. Finally, super- 
shear Rayleigh waves may very well propagate at the surface 
of soft media. 

Materials and methods. A plastic sphere of diameter 
6 mm and density 1.2 g/cw? impacts normally the surface 
of a liquid foam, whose instantaneous vertical deformations 
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Figure 1. Map of the foam's vertical displacement, (a) Sketch of the 
experimental setup (0.8 < V < 6.5 m/s, d = 6 mm, H = 65 mm and 
D= \ \5 mm), (b) Example of deformed fringe pattern (correspond- 
ing height map shown in (d)). Panels (c), (d) and (e) show foam sur- 
face height maps after impact for V = 1.82 m/s at times t = 2.3, 8.9 
and 23.5 ms, respectively. 



are measured by means of a space-time resolved profilom- 
etry technique, as shown in Fig. [TJa) and (b). We use 
the well-characterized "Foamy Regular" cream from Gillette 
lfT4irT5l[T9ll2Qll . The dimensions of the cylindrical container 
avoid side effects: its diameter Z) = 11.5 cm is much larger 
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(a) y < y* (b) y > v* 

Figure 2. Spatiotemporal analysis of surface height in two cases: (a) V < V* (V = 1.82 m/s). (b) V > V* (V = 5.53 m/s). Top: spatiotemporal 
diagrams of h(r,t), with time t in horizontal direction and distance to impact center r in vertical direction. Bottom: h(rj) at various distances r, 
gradually increasing from one sphere diameter d (red curve) to 4d (blue curve). 



than the typical observed attenuation length A (about 1 cm - 
see Fig. [3] (b)) and its depth // = 6.5 cm is larger than A as 
soon as / > 30 Hz ('deep foam' region). After filling the 
tank, a waiting time of 30 minutes is observed for the foam 
to reach a stable state of density 77 g/cm^ . The profilometry 
apparatus consists in a high resolution video projector (Ep- 
son TW5500), projecting a one dimensional fringe pattern on 
the foam surface, and a fast camera (Phantom v9.0) record- 
ing its deformations at 5000 fps. Using the parallel optical 
configuration as in Fourier Transform Profilometry technique, 
the displacement map can be deduced from the local phase 
changes of the fringe pattern [21 1 knowing the optical axes 
separation (23.0 cm), the entrance pupil height (98.5 cm) and 
the projected fringe wavelength (1.82 mm). For better spatial 
accuracy, the phase extraction algorithm used Empirical Mode 
Decomposition filtering 1221 instead of the usual Fourier fil- 
ter |21 1. The typical height resolution is 75 /im (0.05 rad 
in phase) at each point of a grid of 225 /im mesh size and 
the temporal resolution is 0.2 /is. Examples of instantaneous 
height maps are shown in Fig. [T] and in supplementary video 
2. 

Results. Figures [TJc-e) display experimental height maps 
of the sample surface i.e., both the foam and the top of the 
sphere, for an impact velocity of V = 1.82 m/s. Shortly af- 
ter impact, a small bump (about 100 fim amplitude) propa- 
gates radially (Fig. [TJc), ^ = 2.3 ms). The sphere entering 
the foam then creates a dynamic crater, visible on Fig. [TJd) 
(r = 8.9 ms), that disappears on Fig. [TJe) (t = 23.5 ms) after 
the sphere "bounces" and its top becomes visible again above 
the foam surface. One can immediately notice that this de- 
formation is axisymmetric around the impact point. We av- 
erage measurements over all orientations to obtain h(r^t) as 
shown in Fig. [2|a) in a spatiotemporal diagram of this very 



experiment. The amplitude of deformations in this particular 
case is about 1 mm and remains well below d in all our ex- 
periments, provided that observation is made far enough from 
the impact point, at a distance r > 0.7 5d. For clarity, in the 
bottom image, we isolate 4 distances (J, 2 J, 3d and 4J), for 
which we plot the time evolution h(t). For r = d, the signal 
h(t) looks like a strongly attenuated oscillatory motion, with 
a period of approximately 30 ms. At larger distances the sig- 
nal has a similar shape, but delayed and damped. These four 
curves comparison provides evidence for damped propagative 
surface waves with a typical wave speed of a few meters per 
second and centimetric attenuation length. At an higher ve- 
locity of 5.53 m/s, the shape of h(t) is qualitatively different, 
as illustrated in Fig. [2jb). The second maximum present in 
Fig.[2ja) disappears and the long time dynamics is then similar 
to that of an underdamped oscillator. We notice that, surpris- 
ingly, while the impact velocity is multiplied by 3 and there- 
fore the kinetic energy by 9, the crater at r = d remains of 
comparable depth. 

Using wave analysis tools, we now characterize each of 
these two regimes. Adjusting the curves of hf(r) and $/(r), 
the amplitude and phase of the temporal Fourier transform of 
/z(r,0, we can measure the phase velocity V0(/) and attenua- 
tion length laif) defined by 

Both V(f)(f) and /«(/) dispersion curves are shown in Fig.|3] 
The complex formalism allows us to simultaneously account 
for the elastic and dissipative effects by defining the complex 
wave number kdf) = Vphi(f)/27Tf-i/la(f). The wavelength 
A = 27r/5R(^c) and attenuation length la are found to be of 
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Figure 3. Complex dispersion relation components, (a) Phase veloc- 
ity, (b) Attenuation length. Black lozenges: slow impact case. Red 
squares: fast impact case. The error bars correspond to statistics over 
12 experiments of distinct impact speeds and foam realizations. Blue 
line: predictions of shear wave speed V5 based on published law for 
Gillette foam's shear modulus 1 191. 



about 1 cm, much larger than the bubble size (22 /im) and 
smaller than the tank dimensions. Therefore, we consider the 
foam as a semi-infinite homogeneous viscoelastic medium, 
whose storage and loss moduli may depend on frequency. At 
the surface of such a material, the velocity of Rayleigh waves 
is a function of Poisson modulus v and bulk shear wave 
velocity V5', and is always smaller than 1231 . In the case 
of an incompressible material such as the foam used in this 
experiment, v = 0.5 and VR(f) - 0.955 Vs(f) |5 1. In the slow 
impact case, increases from 2 m/s at 40 Hz to 4 m/s at 
300 Hz. The shape of V^C/) curve is also found independent 
on the sphere's diameter. The blue line in Fig. [3] represents 
predictions of the velocity Vs(f) of shear waves based on a 
published law for Gillette foam's shear modulus 1191 . The 
fact that V(f)(f) is very close to Vs(f) strongly suggests that 
the observed waves are Rayleigh waves. Dissipative counter- 
parts of the Rayleigh waves propagating at solid elastic sur- 
faces are called quasi-elastic Rayleigh waves f23\ . In the fast 
impact case, we detect a significantly higher phase velocity 
when / > 190 Hz. Surprisingly, some waves travel much 
faster than bulk shear waves. The reason of this increase is 
ambiguous in the framework of monomodal Fourier analysis: 
either a non-linear increase of elastic modulus modifies the 
speed of shear and Rayleigh waves, or the superposition of 



a second, faster, wave creates a larger apparent velocity. By 
further increasing impact velocity, we obtain the h(t) curves 
shown in Fig. |4] (a) where two bumps propagating at differ- 
ent velocities are clearly distinguishable. This allows us to 
eliminate the first option. The slower bump speed is compara- 
ble to the quasi-elastic Rayleigh wave speed, while the other 
bump travels roughly twice as fast. For r < 3d, the fast bump 
predominates, and the monomodal analysis can be applied for 
measuring its dispersion relation (Fig.|4](b)). If / < 200 Hz, 
it again coincides with the quasi-elastic Rayleigh wave dis- 
persion relation but for larger frequencies jumps to nearly 
twice the expected velocity. 

Discussion. We now discuss the nature of these fast surface 
waves. A liquid-like wave is first considered. Soft interfaces 
can mix elastic and liquid surface properties as in agarose gel, 
on which a pseudo-elastic wave and a pseudo-capillary wave 
can coexist 1241 . For comparing the elastic-like and liquid- 
like contribution we define the ratio of the surface stress in- 
duced by elastic surface waves and liquid surface waves 1 25| : 
P = i^k+ pg/k)/ ji with (1 the shear modulus and 7 the surface 
tension. For the agarose experiments |24l [25l /3 > 1 while 
for liquid foam we find P < 0.02. We conclude that liquid- 
like surface waves do not match our observations. A sec- 
ond possible nature is a different kind of elastic surface wave: 
Rayleigh waves are the real solution of a secular equation, for 
which a complex solution has also been proposed, for elas- 
tic ('leaky waves') (SlIHl or viscoelastic ('viscoelastic waves') 
half- spaces |TOl [TT] |23l. For unity, we name complex solu- 
tions 'supershear Rayleigh waves'. They correspond to waves 
propagating faster than the shear waves and weakly transverse 
in polarization. They may exist only in a limited range of 
Poisson modulus u and shear wave quality factor, defined as 
Q = -^(k^,)/^(k^)\ V > 0.26 in any case and Q < 6.29 

for the viscoelastic case |23 1. However, the physical admissi- 
bility of this solution has been debated |7][TT||261. It has been 
rigorously proven spurious for plane waves but a careful ex- 
amination of the radiation condition for a source localized in 
time and space (as for geological events or impacts) showed 
that it can then be acceptable |9|, and a few fast waves propa- 
gating on soils have also been reported | 27, 28 1. With v = 0.5 
and 2 1 .32, our experiments are well inside the theoretically 
expected existence domain of supershear Rayleigh waves. To 
go further, we predict the supershear Rayleigh wave velocity, 
starting with the measurement of Vr in the slow impact case. 
If VsR is the complex wave speed of the supershear Rayleigh 
wave, then for an incompressible medium, 

(3 .54 + 2.230 V| L23J • Eliminating Vs leads to 



0.91 y| and 



: (3. 88 + 2.44/) Vi- 
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Measurements of Vr and subsequent predictions of Vsr are 
plotted in Fig.|4jb) (blue and red continuous lines). While the 
phase velocity after the 18.3 m/s impact is still close to Vr for 
/ < 200 Hz, it jumps to values in good agreement with the su- 
pershear wave speed Vsr for higher frequencies. This is strong 
evidence that we observed the propagation of the supershear 
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Figure 4. Evidence of two different waves after a very fast impact 
(5 mm diameter bead, V = 18.3 m/s). (a) Detail of h(r,t) at early 
times for values of r ranging from 2d to 4d. (b) Black squares: phase 
velocity dispersion relation (Error bars corresponding to the confi- 
dence interval of the linear regression of (/)(r)). Blue continuous line: 
dispersion relation of the anelastic Rayleigh waves. Red continuous 
line: Prediction of the supershear wave velocity from the slow impact 
dispersion relation, using [T] The confidence interval (dashed lines) 
is defined as in Fig.|3] 



Rayleigh waves on a foam surface. They were visible for im- 
pact velocities faster than the measured shear wave velocities 
and only for / > 200 Hz. This might be due to a deep tank 
condition (experimentally A/if<0.57±0.07) related to the 
need of a companion bulk shear wave | 9 |. 

We conclude that supershear Rayleigh waves have been ob- 
served in our experiment. This implies that they might also 
be expected to propagate on other materials of similar Q. In 
particular, some water or oil filled soils | 29 1 of high seismic 
hazard 1301 , or frictionally held soils. The existence of supers- 
hear waves had also been predicted, in the purely elastic case, 
i.e. without any restriction of Q [9 |. If so, this would imply 
that these waves may be relevant for any type of high Poisson 
ratio soil, not only in soft dissipative ones. 

The authors wish to thank Agnes Maurel, Philippe Petit- 
jeans and Vincent Pagneux for giving them access to the pro- 
filometry setup, Andrea Laugere and Babacar Gomis for their 



contribution to the preliminary experiments, Marie Le Mer- 
rer, Pablo Mininni and Meredith Root-Bernstein for careful 
reading of the manuscript, as well as Etienne Couturier, Paul 
Martin and members of the GDR Mousses for inspiring dis- 
cussions. 



[1] M. Bouchon and M. Vallee, Science 301, 824 (2003). 
[2] K. Xia, A. J. Rosakis, and H. Kanamori, . 
[3] L. Freund, Dynamic Fracture Mechanics (1998). 
[4] S. Das, Science 317, 905 (2007). 
[5] L. Rayleigh, R Lond. Math. Soc. sl-17, 4 (1885). 
[6] R Gilbert and S. J. Laster, B. Seismol. Soc. Am. 52, 299 (1962). 
[7] M. Hayes and R. S. Rivlin, Z. Angew. Math. Phys. 13, 80 
(1962). 

[8] J. A. Hudson, The Excitation and Propagation of Elastic Waves 
(1980). 

[9] C. T. Schroder and W. R. S. Jr., J. Acoust. Soc. Am. 110, 2867 

(2001) . 

[10] J. G. Scholte, Physica 13, 245 (1947). 

[11] P. K. Currie, M. A. Hayes, and P. O'Leary, Quart. Appl. Math. 
35 (1977). 

[12] M. Romeo, Appl. Math. Lett. 15, 649 (2002). 
[13] R. Hohler and S. Cohen-Addad, J. Phys.: Condens. Matter 17, 
R1041 (2005). 

[14] N. Mujica and S. Fauve, Phys. Rev. E 66, 021404 (2002). 
[15] A. D. Gopal and D. J. Durian, J. Colloid Interface Sci. 213, 169 
(1999). 

[16] P. J. Cobelli, A. Maurel, V. Pagneux, and P. Petitjeans, Experi- 
ments in Fluids 46, 1037 (2009). 

[17] A. Maurel, P. Cobelli, V. Pagneux, and P. Petitjeans, Applied 
optics 48, 380 (2009). 

[18] G. Lagubeau, M. Fontelos, C. Josserand, A. Maurel, V. Pag- 
neux, and P. Petitjeans, Physical review letters 105, 184503 
(2010). 

[19] K. Krishan, A. Helal, R. Hohler, and S. Cohen-Addad, Phys. 

Rev. E 82, 011405 (2010). 
[20] D. J. Durian, D. A. Weitz, and D. J. Pine, Phys. Rev. A 44, 

R7902 (1991). 

[21] M. Takeda and K. Mutoh, Appl. Opt. 22, 3977 (1983). 

[22] N. E. Huang, Z. Shen, S. R. Long, M. C. Wu, H. H. Shih, 

Q. Zheng, N.-C. Yen, C. Tung, and H. H. Liu, Proc. R. Soc. 

A 454, 903 (1998). 
[23] J. M. Carcione, Wave Fields in Real Media: Wave Propagation 

in Anisotropic, Anelastic, Porous and Electromagnetic Media, 

Vol. 38 (2001). 

[24] F Monroy and D. Langevin, Phys. Rev. Lett. 81, 3167 (1998). 
[25] K. Ahn, Y. Kosevich, and M. W. Kim, Europhys. Lett. 60, 241 

(2002) . 

[26] C. G. Lai and G. J. Rix, B. Seismol. Soc. Am. 92, 2297 (2002). 
[27] M. Roth, R. Spitzer, and F. O. Nitsche, Proceedings of 12th 

SAGEEP, Oakland, U.S.A. PROCEEDING , 49 (1999). 
[28] J. S. Martin, G. D. Larson, and W. R. S. Jr., J. Acoust. Soc. Am. 

120, 2676 (2006). 
[29] V. A. Korneev, G. M. Goloshubin, T. M. Daley, and D. B. Silin, 

Geophysics 69, 522 (2004). 
[30] H. Flores-Estrella, S. Yussim, and C. Lomnitz, Nat. Hazards 

40, 357 (2007). 



